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INVESTIGATION OF SOME TURBUIENT-BOUNDAEY -LAYER 


VELOCITY PROFILES AT A TUNNEL NALL 


WITH MACH NUMBERS UP TO 1.2 


By Marshall P. Tulin and Ray H. Wright 


SUMMARY 


Turbulent— boundary— layer profiles at large Reynolds numbers and 
with Mach numbers up to 1.2 are presented. Under the conditions of this 
investigation, the velocity profile shapes were substantially unaffected 
by compressibility. The nondimens ional profile shape was found to be a 


displacement thickness and momentum thickness, respectively, which would 
be obtained if the velocity profiles were analyzed as for the incompressible 
f lowj . At a Mach number of 1.2 and with the low values of H^ existing 

for this investigation, the turbulent boundary layer under the shock did 
not separate. The increase in displacement thickness across the shock 
could be calculated by use of the momentum equation. 


An investigation of turbulentHDOundary— layer profiles at large 
Reynolds numbers and Mach numbers up to 1.2 was carried out in connection 
with the design and investigation of a circular supersonic nozzle for 
Mach number 1.2. No systematic variation of the variables influencing 
the boundary— layer development was attempted. Only the boundary— layer 
profiles existing at various positions on the nozzle wall were investigated 

The results obtained in this investigation are turbulent— boundary- 
layer— profile data at large Reynolds numbers, of the order of 40,000 
based on the momentum thickness, and at Mach numbers up to 1.2. The 
effects of compressibility on the velocity profile shapes and form 
parameters were investigated. Some boundary— layer shock— interaction 
data are presented. 

The experimental data obtained are of considerable value, particu- 
larly because no experimental work investigating the effect of compressi- 
bility on the turbulent— velocity profile Bhapes has previously been 


function only of the parameter H^ 



where 5^ and 9 ^ are the 
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reported, although semiempirical methods for calculating the thickness 
of the incompressible turbulent boundary layer hare, by the use of 
reasonable assumptions, been extended to compressible flow (reference 1) . 


STMBOIS 


H 

c 


H 


i 


M 


E 0 


c 


T o 



ratio of displacement thickness to momentum thickness for 
compressible flow 

ratio of displacement thickness to momentum thickness for 
incompressible flow 

Mach number outside boundary layer 

Reynolds number based on boundary— layer momentum thickness 



absolute stagnation temperature inside boundary layer 
absolute stagnation temperature at outer edge of boundary layer 




u 

U 

x 

7 

5 

V 


velocity inside boundary layer 
velocity outside boundary layer 
distance from throat, positive downstream 
distance normal to surface 

distance from surface to outer limit of boundary layer 
boundary— layer displacement thickness for compressible flow 
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6 i* 
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0 i 


P 


■boundary— layer displacement thickness for incompressible flow 

T/o ( i_ ^ 

boundary— layer momentum thickness for compressible flow 

(/o 

boundary— layer momentum thickness for incompressible flow 

(f>^) 

coefficient of viscosity at outer edge of boundary layer 

density inside boundary layer 
density at outer edge of boundary layer 


APPARATUS AND METHODS 


This investigation was carried out in the Langley 8— foot high-speed 

tunnel. The Reynolds number per foot was of the order of 3 x 10 ^ and 
varied from l8,000 to 64,000 based on the boundary— layer momentum thick- 
ness. The variation in Mach number was continuous up to nearly 1.00, 
but only one supersonic Mach number, approximately 1.2, was attainable 
with approximately zero pressure gradient at the measuring station. 

Mach numbers between 1.00 and 1.2 existed only in the accelerating flow 
between the throat and the supersonic test section. Except for the 
possibility of moving to various positions' along the nozzle, the 
boundary— layer Reynolds number was not varied independently of the Mach 
number, and only two regions existed, one near the throat and the other 
at the supersonic test section, where the pressure gradients were 
approximately zero. 

The principal measuring stations are indicated in figure 1, which 
also shows the nozzle profile shape along with representative Mach 
number distributions. The data presented in figures 2 and 3(a) were 
obtained at these stations. In these figures identification of the 
test conditions is facilitated by the values of x and M given. It 
will be noted that for the data of figures 2 and 3(a) long regions of 
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accelerating flow preceded the test, positions. For the general analysis 
of figure 4, data from various positions on the wall of the nozzle and 
also from a previously used subsonic test section are included. The 
data presented in figures 5 and. 6 were obtained under the normal shock 
terminating the supersonic flow at the downstream end of the supersonic 
test section (see fig. l) . 

In order to find out whether the boundary— layer development was 
affected by condensation of water vapor, measurements were taken with 
the tunnel sufficiently cool to cause considerable fog in the flow, and 
these measurements were compared with those taken at much higher 
temperatures. No appreciable difference appeared in the boundary— layer 
measurements . 

The methods of reference 1, in combination with some boundary- 
layer measurements in the entrance cone of the tunnel, were employed 
in computing the boundary— layer displacement thickness used in the design 
of the supersonic nozzle. The success of the nozzle designed in 
producing the design supersonic flow attested to the over-all accuracy 
of the boundary— layer calculations; but because of local pressure 
gradients at the walls and because of asymmetry angularly about the 
tunnel wall of the wall boundary layer (the cause of this asymmetry has 
not been determined) the boundary— layer calculations could not be 
accurately checked by means of boundary— layer surveys along a single 
axial line. The present paper is therefore limited to the investigation 
of boundary— layer profiles and of the behavior of the turbulent boundary 
layer -under the terminal normal shock. 

The total pressures through the boundary layer were measured by 
means of banks of tubes open to the oncoming flow. The banks were of 
two sizes extending outward into the flow 3 inches and 6 inches, respec- 
tively. For the first 3 inches the total-preBsure tubes were spaced 

at inch intervals measured from center to center. On the 6— inch 
^ 1 

banks the tubes were spaced at — —inch intervals, measured from center to 

center, over the outer 3 inches. The pressure differences were measured 
and photographically recorded on manometers containing tetrabr omoe thane . 
Only data near a Mach number of 0.5 and above are presented. The total- 
pressure tubes used were of 0-.050— inch— outside-diameter tubing. In a 
few cases the total— pressure tubes were interspersed with total 
(stagnation) temperature probes constructed by placing small thermocouple 
heads at the entrances to total— pressure tubes; 0.060— inch— outside- 
diameter tubing was used for these probes. Temperature probes calibrated 
in the free stream where the stagnation temperature is known showed only 
small corrections required. With the small supersonic Mach numbers 
encountered in this investigation, the corrections to the total pressure 
were small. 
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For taking accurate measurements of the boundary— layer profile at 
a given point oyer a period of time, single total— pres sure , static- 
pressure , and temperature probes were mounted on a micrometer screw. 

With this apparatus the position of the tubes was continuously variable 
from the wall outward .to 3 inches from the wall. 

For most of the surveys the static pressure was obtained from 
measurements of the pressure at orifices in the tunnel wall. These 
orifices were 0.031 inch in diameter. Except in the region under the 
terminal normal shocks the pressure gradients were sufficiently small 
to cause the static pressure to be essentially constant through the 
boundary layer. 

Mach numbers were obtained from the ratio of static pressure to 
total pressure. Yelocity and density were computed under the assumption 
of constant stagnation temperature through the boundary layer. Under 
operating conditions the stagnation temperature in the tunnel normally 
rises as much as lkO° F, but because of inadequate mixing of the cooling 
air, the temperature near the wall was as much as 50 ° less than that 
near the center of the tunnel. 

The boundary— layer data have been reduced as if they had been 
obtained on a flat rather than a cylindrical surface. Because the 4— foot 
radius is sufficiently large, the error due to this simplication in 
calculation of the momentum thickness 6 C and of the displacement 
thickness 5 C * is less than 2 percent and that in the form param- 

eter H = is negligible. 

o y c 


EESULTS AND DISCUSSION 


Temperature considerations .— Temperature data obtained were 

incomplete and not sufficiently reliable to permit the exact calculation 
of velocities within the boundary layer and of the boundary-: layer 
displacement and momentum thicknesses. . The values of these quantities 
presented herein are therefore based on the assumption that the stag- 
nation temperature is constant throughout the boundary layer. Some 
stagnation— temperature profiles, which, however, are not regarded as 
completely reliable because each profile was taken point by point over 
a period of time during which steady temperature conditions may not 
have existed, indicated a maximum change in absolute stagnation temper- 
ature throughout the boundary layer not exceeding ±6 percent. At a 
Mach number of 1.2, the largest for these tests, the absolute stagnation 
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temperature exceeds the stream temperature by somewhat less than 30 per- 
cent of the stream temperature; and if, as has been found by other inves- 
tigations, with a turbulent boundary layer approximately 90 percent of 
this temperature difference is recovered at the wall provided no heat 
transfer takes place through the wall, the change in stagnation temper- 
ature through the boundary layer at this Mach number is about 3 percent. 
With these considerations in mind, a linear variation of stagnation 
temperature amounting to a 6-percent decrease at the wall was assumed, 
and the resulting errors involved in assuming constant stagnation tem- 
perature were estimated. The following percentage errors were found in 
the various quantities: 


Maximum error 
(percent) 

.... 3 

• • • • — 

.... —7 

.... —2 

.... 7 

.... —10 

.... 19 

The probable errors are, of course, much less than these maximum values. 
In particular the values of u and presented herein are not believed 

to be appreciably in error on account of the stagnation— temperature 
variation. 

Telocity profiles .— Typical velocity profiles with Mach numbers 
from 0.5 to 1.19 are shown in figure 2. The velocity within the boundary 
layer u in terms of the velocity U outside the boundary layer is 
plotted against the nondimens ional distance normal to the wall jJ9j> 

where 0 j is the value of the boundary— layer momentum thickness that 
would be obtained from a plot of u/U against y if the flow were 
assumed inconpressible . 

With incompressible flows, all nondimensional turbulent— b oundary— 
layer profile shapes have been found to be approximately functions of 
the single parameter 


Quantity 

u • • 

8 i ’ • 

0 i * • 
H ± . . 

8 c * * 
0 c * * 
H c * * 




( 1 ) 
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(See, for instance, reference 2.) It is of interest to discover 
whether a similar relation holds also with compressible flow. The com- 
pressible nondimens ional velocity profiles have therefore been analyzed 
as .if they applied to incompressible flow. The values of 0^, and 

their ratio 1L have been determined, and values of u/U at various 
values of y/0^ have been plotted against in figure 4. The symbols 

apply to the data of this investigation, which include measurements at 
Mach numbers from 0.5 to 1.2 with variations in Reynolds number caused 
by changes in speed and variations in position along the tunnel wall. 

In some cases, the pressure gradients immediately upstream from the 
measuring position are positive and in some cases negative. Measurements 
in the boundary layer under normal shock are also included. The solid 
lines in figure 4 are taken from the low-speed data of reference 2. 

Since all the u/U values at each value of y jQ^ arrange themselves . 

approximately into single lines, that is, a definite value of u/U 
corresponds to each value of y/fL for each value, it may be 

concluded that the nondimens ional velocity profile shape in compressible, 
as well as in incompressible flow, is approximately a function of the 
single parameter H^. Furthermore, because the data of this paper agree 
approximately with the low-speed data of reference 2 (note approximate 
agreement between lines and symbols in fig. 4), the functional relation 
between and the profile shape is the same for compressible as for 


incompressible flow, and thus could be approximately represented by a 

u /y\l/m 

formula of the type - = ! M as in reference 1. 


Most of the values of BL found in this investigation are smaller 

than those reported in reference 2. (See fig. 4.) The reason for the 
small values of obtained is believed to be the existence of the 

-long region of accelerating flow which, in most cases, preceded the 
measuring point. The large Reynolds number R@ is also favorable to 

low values of . Values of EL as great as 1.3 were obtained only 
in the region of shock. 

Because the profile shapes depend on H^, the effects of compressi- 
bility on the profile shapes can be examined by investigating the varia- 
tion of that parameter with Mach number. No exact study of the variation 
of with Mach number can be made from the data of this investigation, 

because the conditions of the tests could not be so controlled as to 
maintain constant other variables, such as pressure gradient and 
boundary— layer Reynolds number, on which the development of might 

depend. The principal survey stations (fig. l) were specially chosen. 
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however, to minimize these extraneous effects, and values of 

obtained from velocity profiles at these stations have been plotted 
against Mach number in figure 3(a). The values of for Mach numbers 

less than 0.9 were obtained at a single station and, in all cases, 
similar flow conditions prevailed. The pressure ' gradients at the survey 
stations were nearly zero and each of the three survey stations was 
preceded by a long region of accelerating flow (see fig. l). The survey 
station for a Mach number of 1.19 was a considerable distance downstream 
of the other two stations and was, therefore, preceded by additional 
accelerating flow. This condition may account for the fact that the IL 

value at that station is somewhat less than the values found at the 
upstream positions. With incompressible flow and with the pressure- 
gradient conditions existing for these data, such relatively small 
values of Hj_ change very slowly and it is therefore reasonable to 
suppose that in the present case these conditions exert only a minor 
influence. It seems unlikely in any case that Hp variations due to 

other causes could approximately compensate any variations due to 
changes in Mach number. The approximate constancy of the Hj_ values 

in figure 3(a) is therefore taken as an indication that and hence 

also the nondimens i anal velocity profile shape are not greatly affected 
by compressibility; that is, for the conditions of the present 
investigation the profile shapes are essentially unaffected by 
compressibility. 


Important because of its occurrence in the momentum equation used 
in the boundary— layer calculations is the quantity 



( 2 ) 


where & c * is the actual boundary— layer displacement thickness and 9 C 
is the actual momentum thickness: 


6 c* 







(3) 

(b) 


Because in these equations the velocity ratios u/U at fixed values 
of y/0 i are unaffected by compressibility, H c must be dependent on 

the density ratios P^Pg. Tn the absence of a pressure gradient through 

the boundary layer, the density depends only on the temperature, which. 
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with an insulated boundary. Increases from stream temperature outside 
the boundary layer to a value almost equal to Btream stagnation tempera- 
ture at the wall. Because the difference between stagnation temperature 
and stream temperature increases with Mach number, the density ratio and 
thereby also H c depend on the Mach number. This dependence is shown 

in figure 3(b) where, for a velocity profile 



that yields a value of IL equal to about 1.25, H c is shown plotted 
against Mach number. For the top curve (fig. 3(h)) the stagnation 
temperature is assumed constant. 

In order to show the sensitivity of H c to the temperature distri- 
bution through the boundary layer, values of H c have been computed for 
two cases of cooling at the wall. In both cases, a temperature at the 
wall 6 percent less than the absolute stagnation temperature of the 
stream has been assumed. Stagnation— temperature distributions through 
the boundary layer have been chosen similar to velocity profiles 

ft m °' 94 + °’ 06 (I) 

5 

Variations of H c for n = 5 and- n = 8 are shown (fig. 3(h)). The 
cooling has a substantial effect on the values of H c , though less than 

if the linear temperature distribution used for estimating the possible 
errors due to neglecting the stagnation— temperature variation had been 
assumed. 

From equations (3) and (4) it is evident, because the density within 
the boundary layer is decreased on account of the temperature increase, 

that with a given velocity profile against jj the displacement 

thickness S * is increased and the momentum thickness 9 n is decreased 
c c 

by compressibility. With a constant value of 0 C , which is the usual 

condition for flat-plate flow (see discussion and references of refer- 
ence 1 concerning absence of compressibility effect on skin friction), 
the displacement thickness must increase still more with increase in 
Mach number. Thus, the usual effect of compressibility on the turbulent 
boundary layer is to increase its displacement. 
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Turbulent boundary layer under normal shock .— Figure 5 shows the 
variation of indicated (from wall pressures) Mach number M, displace- 
ment thickness & c *, and form parameter through the region under 

a normal, shock near the end of the supersonic test section. In order 
to facilitate the experiment , the shock wave was moved past the bank of 
survey tubes rather than vice versa. This process renders the Mach 
number inaccurate, because for downstream points in figure 5 the Mach 
number preceding the shock is actually lower than that shown. Never- 
theless the phenomena shown are believed to be qualitatively correct, 
and the experimental method used has the advantage that the effects 
shown occur at a given station and variations due to movement from one 
position to another in the tunnel are eliminated. Because the only 
changes involved are changes in Mach number distribution, the variations 
shown are due solely to primary or secondary effects of changes in Mach 
number and Mach number gradient. Because of the very thick boundary 
layer (& c * £ 0-5 inch) and, also perhaps, partly because of the observed 
unsteadiness in shock position, the region of large pressure gradient at 
the wall is spread over a distance of some 15 inches. 

The sharp increase in displacement thickness with decrease in 
velocity as the shock Is moved upstream past the survey position should 
be noted (fig. 5)-' This thickening is, of course, to be expected. By 
assuming that the inertia terms were predominant for the flow through 
the shock and by choosing from figure 6 a value of H c , which in this 

case was constant through the shock, the momentum equation could be 
used for computing the boundary— layer displacement thickness. True, 
some of the approximations involved In that equation, particularly as 
to constancy of pressure through the boundary layer and thinness of the 
region of pressure jump under the shock, are not well satisfied. Never- 
theless, a point so computed (approximate theory in fig. 5) falls near 
the experimental curve. It seems quite likely, therefore, that with 
Reynolds numbers sufficiently large to insure turbulent boundary layers 
the approximate theory may serve as a useful guide to the behavior of 
the boundary layer under shock. 

As may be observed from figure 5* "the quantity H^, which charac- 
terizes the shape of the velocity profiles, undergoes the expected 
increase in the region of positive pressure gradient under shock. The 
variation in the parameter H accompanying this increase in IL 
through the shock' should not be expected to be as great because of the 
effect of the decreasing Mach number In reducing H c . The exact 

behavior of H c is dependent on the stagnation— temperature profile. 
Figure 6 shows the variation of values of H c calculated on the basis 

of a constant stagnation temperature throughout the boundary layer. 

From this condition, H c is observed to be nearly constant throughout 
the region of shock. 
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Because of the relatively small values of the supersonic Mach 
number attained, the variations of the houndary— layer characteristics 
in the region under the normal shock were less than those reported in 
reference 3, for which Mach numbers up to 1.4 existed. In the present 
investigation no evidence of separation was found, though special- 
efforts were made to detect it. With higher Mach numbers and, con- 
sequently, more severe shocks, separation is to be expected at least 
locally . 


CONCLUSIONS 


From the results of an investigation of turbulent— boundary— layer 
profiles at large Eeynolds numbers and with Mach numbers up to 1.2, the 
following conclusions are made: 



1. At low values of the velocity-profile— shape parameter H. 

where Eh and 0 i are the displacement thickness and momentum 

thickness, respectively, which would be obtained if the velocity profiles 
were analyzed as for incompressible flow'j found in this investigation, 

no important effects of compressibility on turbulent— boundary— layer 
profile shapes at large values of Eeynolds numbers were found. 


2. With compressible flow, as with incompressible flow, the profile 
shape is a function of the single parameter 


3« With a given value of the ratio H c of actual displacement 

thickness to actual momentum thickness is a function of Mach number and 
of the stagnation— temperature distribution through the boundary layer. 


4. Because of the heating within the boundary layer, its displacement 
thickness is increased with increase in Mach number. 


5 . The pressure jump through the shock is spread over a distance 
at the surface that is large compared with the boundary-layer displace- 
ment thickness. 


6. With the comparatively weak normal shock corresponding to a 
Mach number of 1.2 and with a comparatively low value of H i upstream 
of the shock, separation of the turbulent boundary layer did not occur. 
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7 . For the conditions of this investigation, the momentum equation 
could he used for calculating the increase in displacement thickness 
across the shock. 


Langley Aeronautical Laboratory 

National Advisory Committee for Aeronautics 
Langley Air Force Base, Ya. 
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Figure 1.— Mach number distribution at nozzle wall. (Symbols represent runs at various speeds. 
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Figure 3 — Variation of the shape parameter H with Mach number. 
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Figure 4*— Comparison of velocity profiles of present investigation with those of reference 2 with 
respect to incompress ihle form parameter 1^. (Symbols apply to present investigation and lines 

to that of reference 2.) 




NACA RM L9H29a 


21 




Figure 5*— Turbulent boundary— layer characteristics In the region of 
shock as measured on the wall of the Mach number 1.2 nozzle in 
the Langley '8— foot high-speed tunnel. Reynolds number per foot 
is 3*8 X 1()6. 
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Figure 6.— Variation of form parameter under shock. 


